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Born and Oppenheimer reported an approximate separation of molecular eigenfunctions 
into electronic, vibrational, and rotational parts, but at the end of their paper showed that 
the two angles describing rotation of the nuclei in a diatomic molecule are exactly separable. 
A year later in a two-part work devoted strictly to diatomic molecules, Wigner and Witmer 
published (1) an exact diatomic eigenfunction and (2) the rules correlating the electronic 
state of a diatomic molecule to the orbital and spin momenta of the separated atoms. The 
second part of the Wigner-Witmer paper became famous for its correlation rules, but, oddly, 
the exact eigenfunction from which their rules were obtained received hardly any attention. 
Using three fundamental symmetries, we give a derivation of the Wigner-Witmer diatomic 
eigenfunction. Applications of our derivations are fundamental to predicting accurate di¬ 
atomic molecular spectra that we compare with recorded spectra for diagnostic purposes, 
such as measurements of molecular spectra following generation of laser-induced plasma. 
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I. INTRODUCTION 


Q 


In the introduction to their paper, Born & Oppenheimer |l( allude to an exact separation of 
two rotational coordinates in the diatomic molecule. In their next section, which is applicable to 
polyatomic molecules, they introduce a coordinate system attached to the nuclei whose orientation 
is set by the Euler angles, and note that there are terms in the molecular Hamiltonian in which 
both electronic and nuclear coordinates appear thereby preventing the exact separation of the total 
eigenfunction into a product of electronic and nuclear eigenfunctions. In their final section, Born 
and Oppenheimer return to the diatomic molecule and give the details of the exact separation of 
two of the Euler angles, angles 9 and <j> that describe rotation of the two nuclei. The spherical 
harmonic Yg m (9,(j)) is the angular momentum part of Born-Oppenheimer diatomic eigenfunction. 
A year after, Wigner & Witmer Q| published a two-part article on diatomic theory in which they 
replace the spherical harmonic Y£ m (0,0) with the Wigner D-function 9, x). The Wiger- 

Witrner paper 2| was the introduction of the rotation matrix elements into diatomic theory. After 
the Wigner-Witmer paper, the Wigner D-function mostly disappeared from diatomic literature 


for about four decades. Hirschfelder Sz Wigner [3j used the ZJ-function in their separation of 6 
coordinates (3 for the total linear momentum, 3 for the total angular momentum) for N particle 
systems, but do not explicitly mention the diatomic molecule. Again not specifically mentioning 

4] repeat the separation of 6 coordinates 


Si 


in his quantum 


the diatomic molecule, Curtis Sz Hirschfelder Sz Adler 

for N particle systems, and consider the three-body system in detail. Davydov 
mechanics textbook used the D-function in his discussion of diatomic theory. At about the same 
time, Rubin 


employed it for his calculations of Honl-London factors. Pack Sz Hirschfelder S 
used the D-function to separate two angular rotational coordinates in the diatomic eigenfunction 
but failed to notice that their Eq. (2.35) holds for all values of the third Euler angle, not just 


n 


7 = 0. Za re | 8p explicitly used the D-function in their case (a) basis function. Judd Q] and 


Mizushima LLOfl, in their treatments of diatomic theory, introduce the D-function and discuss its 
mathematical properties, but do not explicitly display it in their Hund’s cases (a) and (b) basis 
functions. The Wigner D-function has since become a vital mathematical component of diatomic 
theory, as comprehensively collated by Varshalovich Sz Moskalev Sz Khersonskii [ll], However, the 
exact separation of 9 x) in the diatomic eigenfunction where the nuclear coordinates 4> and 

9 appear with the electronic coordinate x has remained essentially forgotten for eight decades. 
We give a derivation of the Wigner-Witmer diatomic eigenfunction based upon three fundamental 
symmetries and the geometrical symmetry of a molecule possessing precisely two nuclei. 
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II. DERIVATION OF THE WIGNER-WITMER DIATOMIC EIGENFUNCTION. 


Here we obtain the Wigner-Witmer eigenfunction by applying three symmetry principles to 
the eigenfunction of a free conservative system composed of N electrons and precisely two nuclei. 
Energy is the generator of translations in time, the time translation (evolution) operator U{t,to ) 
is a continuous unitary operator, the total energy is a constant of the motion, and the dependence 
of the eigenfunction on the physical variable time t is exactly separable if the time origin to can 
be associated with some physical event. Linear momentum is the generator of translations in 
space, the spatial translation operator T(R,Ro) is a continuous unitary operator, the total linear 
momentum is a constant of the motion, and the total linear momentum is exactly separable if 
the coordinates Rcm of the center of mass can be introduced as physical variables of the system. 
Angular momentum is the generator of rotations, the rotation operator 1Z(a, /3, 7 ) is a continuous 
unitary operator, the total angular momentum is a constant of the motion, but the total angular 
momentum 3{(f), 9, x) is no L in general, exactly separable because except for very simple systems 
one cannot find physical rotations cf>, 9 , and x which duplicate the angles a, /3, and 7 of coordi¬ 
nate rotation. The diatomic molecule deserves a special place in the quantum theory of angular 
momentum because it is one of the most complicated systems for which the Euler angles a, f3, 
and 7 of coordinate rotation are also the angles of physical rotation describing the total angular 
momentum J 

Quantum mechanical descriptions of the diatomic molecule typically begin with the Hamilto¬ 
nian, but minutia in the Hamiltonian tend to obscure the few fundamentals at play. For example, 
Brown &: Carrington G write a diatomic Ham ltonian, their Eq. (2.297), containing 32 types of 
Hamiltonian terms. We begin our discussion of diatomic theory with the eigenfunction 


R- 2 , ■ ■ • , Ra; R{ ,,t) — (Rl, R 2 , ■ ■ • ; RjV, Ra; R{» t \nvJM) 


(1) 


in which Ri, R 2 ,..., Rat are the spatial coordinates of the N electrons and R a and Rf, are those of 
the nuclei. The total angular momentum quantum numbers J and M refer to the true total. That 
is, in spectroscopic nomenclature they would be replaced by F and Mp. The symbol n represents 
all other required quantum numbers and continuous indices except the vibrational quantum number 

v. 

The symmetries of translation in time and translation in space produce a separation of the time 
coordinate t and the spatial coordinates Rcm of the center of mass. A two-body reduction of the 
motion of the nuclei requires placement of the coordinate origin a the center of mass of the nuclei, 
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and then replaces R a and Rb with the internuclear vector r. Of the 3N + 7 dynamical variables in 
the total eigenfunction (JTJ), 3N + 3 remain in the internal eigenfunction (ri, r 2 , • • •, r^r, r | nvJM). 
The axes of the translated coordinates ri, r 2 ,..., rjy, r are parallel to those of the original laboratory 
coordinates Ri, R 2 ,..., Rat, R a , Rb- We now address how rotational symmetry influences the 
internal eigenfunction. 

Operation of the rotation operator !Z(a,f3,'y) on the internal eigenfunction yields 

( r i, r 2 ,..., rjv, r \R,{a, (3 : ^)\nvJM) = {r[,r' 2 ,... ,r' N ,r'\nvJM) (2) 


where primes denote rotated coordinates given by 


£>(<*,/3,7) 


cos a cos f3 cos 7 — sin a sin 7 sin a cos f3 cos 7 + cos a sin 7 — sin (3 cos 7 

— cos a cos (3 sin 7 — sin a cos 7 — sin a cos (3 sin 7 + cos a cos 7 sin f3 sin 7 
cos a sin (3 sin a sin /3 cos /3 


( 3 ) 


x' 


X 

y' 

e 

11 

y 

z' 


z 


The effect of lZ(a, (3, 7 ) on the eigenfunction can be rewritten as 


( 4 ) 


(ri,r 2 , • • • ,rjv,r | nvJM) = (r^r^,... ,r' N ,r \Tl\a,/3,^)\nvJM) (5) 

J 

= Y < r 'i> r 2 ’ • • • > (jn \n\a, /3, 7 ) \JM) ( 6 ) 

fl =~J 
J 

= Y ^ r i’ r 2’ ■ • • > r JV> r ' \nvJQ) D-£ n (a,P,j). (7) 

n=-j 

When the spherical coordinates of the internuclear vector r = r(r, 9, 0) are introduced, the equation 
becomes 

J 

(ri, r 2 ,..., r N , r, 6, 4> \nv JM) = Y ( r 'n r 2 , ■ • •, r N, r,Q', <\> \nvJil) D^ n (a, /3, 7 ). ( 8 ) 

n=-j 

The internuclear distance r is unprimed on the right because it is a scalar. Because physical 
rotation 4> and coordinate rotation a are both counterclockwise rotations about the z axis, the 
physical angles q 7 is given by 


( 9 ) 
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Similarly, physical rotation 0 and coordinate rotation (3 are counterclockwise rotations about the 
first intermediate y axis of the total coordinate rotation. 

6* = 6-p. (10) 

Rotational symmetry gives us the option to view the molecule at any orientation we choose, and 
we choose a = 4> and f3 = 6. 

J 

(ri, r 2 , ... ,r N ,r,d,(j) \nvJM) = E ( r i, r' 2 ,..., r' N , r \nvJQ) 9, 7 ). (11) 

n=-j 

The rotated coordinates of the one of the electrons, we arbitrarily select the electron labeled 1, are 
expressed in cylindrical coordinates p \, £{, and x'\ ■ 

J 

(iT, r 2 ,... ,r N ,r,0,<j) \nvJM) = ^ (pi, ([, xi, r 2 , • • •, r 'iv> r \nvjn) D^ Q ((j), 9, 7 ). (12) 

n=-j 

The chosen electron is distance p\ from the internuclear vector and signed distance £[ above 
or below the plane perpendicular to the internuclear vector and passing through the coordinate 
origin. The angle x' describes rotation of this electron about the internuclear distance. Like 
the internuclear distance r, primes on p\ and are unnecessary because they are scalars whose 
values are unchanged by coordinate rotation. Because x\ and 7 are rotations about the same axis, 
coordinate rotation changes, of course, the value of x ’\; but this also means that this angle has a 
value xi br laboratory coordinates. The coordinate rotation angle 7 is chosen to make x\ zero, 

Xi = Xi ~ 7- (13) 

The equation for the eigenfunction now reads 

J 

(p,C,X, r 2 , • ■ ■ ,T N ,r,6,(j}\nvJM) = ^ (p, (, r 2 ,..., rjv, r \nv) D J Mn {^), 9, x) (14) 

n=-j 

after the subscripts on p\, Ci, and xi have been dropped. This result is the Wigner-Witmer 
diatomic eigenfunction. The total diatomic eigenfunction is given as the sum of 2 J + 1 products 
of electronic-vibrational basis functions (p, £, r 2 ,..., r' N , r \nv) and total angular momentum basis 
functions D'^ n ((j),9,x) or Wigner D-functions. The Born-Oppenheimer approximation separates 
the electronic-vibrational basis into the product of electronic and vibrational basis functions, but 
the separation of Df^(<p,9,x) from the electronic-vibrational basis is exact. Many individual 
orbital and spin momenta are contained in the electronic-vibrational basis, but 9, x) is the 

total angular momentum basis. 
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It is noteworthy that eigenfunctions for a rotational state of the diatomic molecule usually 
contain parity, in other words, parity symmetrization is customary. Inclusion of the discrete parity 
symmetry is accomplished after the construction of the Wigner-Witmer diatomic eigenfunction. 
The parity operator can be constructed from proper and improper rotations that have a determinant 
of the transformation matrix of +1 and -1, respectively. Subsequently, Equation 2.14 for the 
eigenfunction can be split for specific values of II followed by parity symmetrization. The approach 
of including parity for specific values of D has been utilized in the literature in order to reduce the 
size of the Hamiltonian matrix prior to finding eigenvalues by diagonalization. In our work |l3]. 
the parity operation is considered after establishment of the eigenfunction in terms of J and M as 
sum over hi in Eq. 2.14. 

Clearly, as electronic states for the diatomic molecule are considered, parity is paramount for 
building these states utilizing the Wigner-Witmer correlation rules [jjJ]. Yet in this work we focus 


on the use of the Wigner-Witmer eigenfunction for computation of spectra, rather than molecular 
structure predictions in non-Born-Oppenheimer calculations for molecules [15, 16]. 


III. DISCUSSION AND CONCLUSION 


The mixing of the electronic coordinate x with the two nuclear coordinates <f> and 9 in 
x) is an obvious departure from current expositions of diatomic theory. However, 
the exact separation of the total angular momentum basis from the electronic-vibrational basis 
proves useful. For example, writing the parity operator V as the product of a proper rotation V a p-y 
and an improper rotation P%, 

V = V a p 7 P e , (15) 

one obtains a simple equation for the parity of diatomic states. The eigenvalue of V is, of course, 
±1, and the product of eigenvalues p-zPap^, 

P = PeR J+2M , (16) 


is always ±1 as required. Sign changes due to parity can show different effects on the Wigner 
D-function. Depending on the specific rotation group, the effect can be expressed in terms of J 
and M or in terms of J and D. Table [I] summarizes the sign changes due to parity. These results 
are consistent with the ones presented by Varshalovich & Moskalev & Khersonskii [111 ]. 

Note that for half-integer J in Equation (3.2), the individual eigenvalues ps and p a p^ are purely 
imaginary. A widely accepted convention allows one to treat the parity eigenvalues ps and p a/ 3 7 as 
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TABLE I. The sign changes on the components x', y' and z' of a coordinate vector r'(x', y', z') produced by 
three different discrete Euler angle transformations, and the effect of these Euler angle transformations on 
D MSii a ,P,i)- _ 

Transformation group Euler angles Coordinates Effect on D^ n (a, 


C 2 (x') 

a —> 7 r + a 

x' 

-)■ x' 



/3 —> 7T — /? 

y' 

-*■ -y' 

{-) J+2U D J M _s l {ot,P,')) 


7 “7 

z' 

-> -z' 


<W) 

a —> 7r + a 

x' 

-»■ —x' 



P 7T — /3 

y' 

-+y' 

(~) J ~ n D^_ n (a, 1 3 , 7 ) 


7 —> 7T — 7 

z' 

-> -z' 


C2W) 

a —> a 

x' 

-»■ —x' 




y' 

-»• -y' 

(-)"° D J M_ n {ot,P, 7 ) 


7 —> 7T + 7 

z' 

->■ z' 



real for both integer and half-integer J [17]]. If one agrees to always subtract 1/2 from half-integer 
values of J, then the diatomic parity can be written as 


P = +PE (-) J 

= -pe (- y - 1 ' 2 


J integer 
J half-integer 


(17a) 

(17b) 


in which ps = ±1 is always real. 

The Wigner-Witmer diatomic eigenfunction simplifies the process in which one infers molecular 
parameters such as the rotational parameter B v and the spin-orbit parameter A v from experimen¬ 
tally measured line positions. Application of our detailed Wigner Witmer eig enfunctions include 


analyses of low- and high-temperature spectra of diatomic carbon spectra 


181 ]. or as another ex¬ 


ample, development of line strengths for specific transitions of the aluminium monoxide (AlO) 
diatomic molecule {3]. The parameters we use are electronic-vibrational matrix elements. In 
turn, the Born-Oppenheimer approximation separates these into the products of electronic matrix 
elements times vibrational matrix elements, and introduces a large number of differential equations 
which couple the many Born-Oppenheimer vibrational states thereby producing a large Hamilto¬ 
nian matrix. Van Vleck transformations reduce the dimension of the Hamiltonian to yield an 
effective Hamiltonian. If in the fitting process, one deals with electronic-vibrational matrix ele¬ 
ments such as B v and A v instead of breaking them into the products of electronic matrix elements 
and vibrational matrix elements, Van Vleck transformations are no longer required. 
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Conversely, the Wigner-Witmer eigenfunction does not reveal how the total angular momen¬ 
tum is built from its components. One must use an angular momentum coupling model which 
has a complete basis. For example, the Hund’s case a basis appropriate to the Wigner-Witmer 
eigenfunction is 

h J + I r * 

|a) = \nJMQ,KSY) = \J - (p, (, r 2 , ■ ■ ■, r ' N ,r \nv) D J Mn {4 >, 0, x) |<ST), (18) 

where = A + S. Written in our notation, current practice replaces the above with 


l«) = £ (p, C, r ' 2 , \nvBoJMnA) (v BO |«>£>£n(6 0, 0) |SE>. (19) 

vbo =0 

The total angular momentum is not exactly separated in this equation and there is the sum over 
a large number of Born-Oppenheimer vibrational states | vbo)- It is often said that the Born- 
Oppenheimer approximation separates the diatomic eigenfunction into electronic, vibrational, and 
rotational states, but this does not hold with spectroscopic accuracy. In the modeling of upper 
and lower Hamiltonians whose term differences accurately agree with measured line positions, one 
must deal with the large set of coupled differential equations that result when Eq. (1191) is inserted 
in the Schrodinger equation. Analytical techniques that employ Van Vleck transformations and 
parity symmetrization have been developed to yield much smaller effective Hamiltonians. 

Except for the simplest of diatomic molecules, ab initio computations are usually not as accu¬ 
rate as results obtained in experimental spectroscopy. For an accurate prediction of a diatomic 
spectrum, one must have values for the molecular parameters such as B v , A v , X v , j v , ..., and their 
centrifugally stretched forms. Computer programs have been developed which find the molecular 
parameters by fitting upper and lower term differences from model Hamiltonians to accurately 
measured line positions. Such a program begins with trial values for the molecular parameters, 
computes theoretical line positions as eigenvalue differences between upper and lower Hamiltonians, 
computes corrections to the trial values of the molecular parameters from the differences between 
the computed and measured line positions, and iterates until corrections to the parameters be¬ 
come negligibly small. There are now many examples of this algorithm for which the errors in the 
computed line positions, i.e., computed vacuum wavenumbers, are not significantly larger then the 
estimated experimental accuracy. With two exceptions, replacement of the Born-Oppenheimer ap¬ 
proximation with the Wigner-Witmer eigenfunction does not significantly alter the flow charts for 
these programs. First, the matrix elements of the effective Hamiltonian are replaced with Hund’s 
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case a matrix elements unmodified by Van Vleck transformations and parity symmetrization. Sec¬ 
ond, the manual enforcement of selection rules is replaced by computation of the Honl-London 
factors {3]. The exact separation of the total angular momentum in the Wigner-Witmer eigen¬ 
function provides simple, accurate computation of the Honl-London factors, and there is but a 
single diatomic selection rule: Transitions for which the Honl-London factor is non-vanishing are 
allowed. Thus, use of the Wigner-Witmer diatomic eigenfunction represents a significant departure 
from current practices in diatomic theory. 
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